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Abstract 

Foliage photosynthetic and structural traits were studied in 15 species with a wide range of foliage anatomies to gain 
insight into the importance of key anatomical traits in the limitation of diffusion of C0 2 from substomatal cavities to 
chloroplasts. The relative importance of different anatomical traits in constraining C0 2 diffusion was evaluated using 
a quantitative model. Mesophyll conductance (g m ) was most strongly correlated with chloroplast exposed surface to 
leaf area ratio (S c /S) and cell wall thickness (T cw ), but, depending on foliage structure, the overall importance of g m in 
constraining photosynthesis and the importance of different anatomical traits in the restriction of C0 2 diffusion varied. 
In species with mesophytic leaves, membrane permeabilities and cytosol and stromal conductance dominated the 
variation in g m . However, in species with sclerophytic leaves, g m was mostly limited by T cw . These results demonstrate 
the major role of anatomy in constraining mesophyll diffusion conductance and, consequently, in determining the 
variability in photosynthetic capacity among species. 

Key words: cell wall thickness, diffusion model, leaf anatomy, leaf structure, photosynthesis, quantitative photosynthetic 
limitations. 



Abbreviations: a, leaf absorptance; (3, fraction of absorbed light that reaches photosystem II; r*, C0 2 compensation point in the absence of mitochondrial respira- 
tion; <f PS », effective quantum efficiency of the PSII photochemistry; AZ_ ias , effective diffusion path length in the gas phase; e PS n, fraction of electrons absorbed by 
PSII; 5, diffusion path tortuosity; /4 mass , photosynthetic capacity per dry mass; A N , net C0 2 assimilation rate; C a , atmospheric C0 2 concentration; C c , chloroplastic 
C0 2 concentration; Q, substomatal C0 2 concentration; C r C c , C0 2 drawdown from intercellular airspace to chloroplasts; D a , diffusion coefficient for C0 2 in the 
gas phase; D L , leaf density; D w , aqueous phase diffusion coefficient for C0 2 ; f ias , volume fraction of intercellular air spaces; F m ', maximum fluorescence in light- 
adapted state; F s , steady-state fluorescence emission; g oe i, partial liquid phase conductance for different portions along cell walls; g cyt , cytosol conductance; g env , 
chloroplast envelope conductance; g ias , intercellular air space conductance to C0 2 (gas phase conductance); g Sq , sum of liquid and lipid phase conductances; 
g m , mesophyll diffusion conductance; g pl , plasma membrane conductance; g s , stomatal conductance to C0 2 ; g to t. total conductance to C0 2 from ambient air to 
chloroplasts; H/(RT k ), dimensionless form of Henry's law constant; J F , linear electron transport rate from chlorophyll fluorescence; J max , maximum photosynthetic 
electron transport rate; K 0 , Michaelis-Menten constant for the carboxylation activity of Rubisco; K 0 , Michaelis-Menten constant for the oxygenation activity of 
Rubisco; / b , biochemical limitation; /_ ch |, length of chloroplasts exposed to intercellular air spaces; L cyt , diffusion pathway length in the cytoplasm; / ias , gas-phase 
limitation; / m , mesophyll limitation; / s , stomatal limitation; M A , leaf mass per area; O, leaf internal oxygen concentration; p„ effective porosity in the given part of 
the diffusion pathway; Q, incident quantum flux density; R, gas constant; fi d , leaf respiration in the dark; r f i , proportional reduction of D w in the cytosol and in the 
stroma compared with free diffusion in water; R L , leaf respiration in the light; S C /o, Rubisco specificity factor; S c /S, chloroplast surface area exposed to intercel- 
lular air spaces per unit of leaf area; S c /S m , ratio of exposed chloroplasts to mesophyll surface areas; S m /S, mesophyll surface area exposed to intercellular air 
spaces per unit of leaf area; S s , cross-sectional area of mesophyll cells in micrograph; SE, standard error; T oh |, chloroplast thickness; T cw , cell wall thickness; 7~ oyt , 
cytoplasm thickness; 7~ k , absolute temperature; 7~ L , leaf thickness; t mes , mesophyll thickness; l/ cm ax, maximum rates for the carboxylation activity of Rubisco; W, 
width of the leaf anatomical section. 
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Introduction 

Leaf anatomical characteristics are key functional and adaptive 
traits determining plant capacity to thrive in specific environ- 
ments, in particular, because these traits also have important 
implications for foliage potential photosynthesis (Niinemets 
et al, 2009a; Scafaro etal, 201 1; Terashima e? al, 201 1). Analysis 
of global variations in leaf functional traits — the leaf economics 
spectrum — has established that the variation in leaf dry mass per 
area (M A ) is strongly associated with other key leaf traits such 
as maximum photosynthetic capacity per dry mass 04 mass ), leaf 
life span, nitrogen and phosphorous contents per dry mass, and 
respiration (Wright et al, 2004). Species with lower M A present 
short leaf life spans, high photosynthetic capacities and nutri- 
ent contents, and low leaf area construction costs, resulting in 
fast growth in environments with high availability of resources. 
In contrast, species with higher M A and lower ^4,^ present the 
opposite suite of traits and have higher cost for leaf area forma- 
tion, particularly due to investment in vasculature and cell walls 
(Niinemets et al, 2007; Hikosaka & Shigeno, 2009) and overall 
improved resistance to low fertility and drought, but low growth 
rates (Niinemets, 2001; Wright et al., 2004). It has been hypoth- 
esized that the negative relationship between M A and photosyn- 
thetic capacity is partly because of greater biomass investment 
in support tissues and cell wall thickening involving stronger 
C0 2 diffusion limitations to photosynthesis (Niinemets, 1999; 
Wright et al., 2004; Niinemets et al, 2007) 

Mesophyll conductance to C0 2 (g m ) is the measure of the C0 2 
diffusion facility from substomatal cavities to the sites of carbox- 
ylation in the chloroplasts (Flexas et al., 2008, 2012) Mesophyll 
conductance is finite and variable and plays a major role in 
constraining photosynthetic productivity (Niinemets et al., 
2009a). Large differences in g m have been shown between and 
within species with different leaf forms and habits (Flexas et al., 
2008; Warren, 2008; Niinemets et al., 2009a, 201 1). Whilst rapid 
changes of g m in response to environmental drivers probably 
depend on biochemical factors such as changes in the permeabil- 
ity of membranes to C0 2 facilitated by cooporins (Hanba et al., 
2004; Flexas et al., 2006, 2012), maximum values of g m for a 
given species or genotype are suggested to be related to leaf ana- 
tomical properties (Niinemets et al., 2009a; Tosens et al., 2012a). 
In particular, it has been shown that leaves with a more robust 
structure and higher M A exhibit lower photosynthetic rates due 
to large C0 2 drawdown from substomatal cavities (Q) to chloro- 
plasts (C c ), Cj-C c , demonstrating that the photosynthetic capac- 
ity is limited by g m (Flexas et al, 2008, Niinemets et al, 2009a). 
Therefore, understanding the structural and physiological basis 
of variation in g m is crucial for understanding photosynthetic 
controls in natural ecosystems and for breeding of plant culti- 
vars with improved photosynthetic characteristics. 

At the leaf level, two components of M A — leaf thickness 
and density — have been proposed to exert opposite effects on 
setting the maximum g m , with increases in thickness increas- 
ing g m and increases in density reducing it (Niinemets et al, 
2009b, Hassiotou et al, 2010). Inside leaves, the C0 2 diffusion 



pathway consists of two phases, an intercellular gas phase and 
a cellular liquid phase, the latter consisting of aqueous and lipid 
components(Niinemets and Reichstein, 2003b; Evans et al, 
2009). The gas phase pathway through intercellular air spaces is 
assumed to have a smaller effect on the overall diffusion limita- 
tions than the components of the liquid phase (Evans et al. , 2009). 
This was confirmed in several studies comparing C0 2 diffusion 
in air and helox — air where helium replaces nitrogen to increase 
diffusivity — showing that the diffusion in the intercellular gas 
phase had little effect on photosynthesis (Parkhurst and Mott, 
1990) The cellular phase is composed of the cell wall, plasma 
membrane, cytosol, and chloroplast envelopes and stroma. 
Among these components, the cell walls and chloroplast enve- 
lope have been suggested to limit g m most severely (Terashima 
et al, 2011). Accordingly, several reports have shown positive 
correlations between g m and the surface of chloroplasts adjacent 
to intercellular air spaces (SJS), which is sometimes considered 
as the most important anatomical trait affecting g m (Evans et al, 
1994; Terashima et al, 2006; Tholen et al, 2008). However, some 
estimates suggest that differences in cell wall thickness (r cw ) can 
explain as much as 25-50% of the variability in g m (Evans et al, 
2009; Terashima et al, 201 1; Tosens et al, 2012b). Negative cor- 
relations between g m and r cw have been shown when comparing 
Australian Banksia species (Hassiotou et al, 2010), rice relatives 
(Scafaro et al, 2011), Eastern Australian species with varying 
anatomy (Tosens et al, 2012b), and Mediterranean Abies spe- 
cies (Peguero-Pina et al, 2012). Recently, Terashima et al (201 1) 
showed that gJ(SJS) decreases with increasing r cw , i.e. the rela- 
tive influence of the exposed chloroplast surface in setting the 
maximum g m is variable, and that this variation can potentially 
be explained by variations in r cw . 

Few previous studies have quantitatively addressed the influ- 
ence of leaf anatomical traits on the diffusion of C0 2 , and these 
studies have focused only on a few species and specific parts of 
the C0 2 diffusion pathway (Evans et al, 1994; Terashima et al, 
2006; Hassiotou et al, 2010; Scafaro et al, 2011; Peguero-Pina 
et al, 2012; Tosens et al, 2012b). Hence, the whole diffusion path- 
way of C0 2 from the substomatal cavities to the chloroplasts has 
not been quantitatively linked to g m in plants with widely varying 
leaf structures and photosynthetic capacities. Furthermore, the 
overall importance of g m in constraining the photosynthetic rate 
in species with different foliage architecture has not been char- 
acterized. To fill this gap, we aimed with the present study: (i) to 
analyse the role of different components of the diffusion pathway 
across a wide range of foliage architectures and leaf photosyn- 
thetic capacities; (ii) to associate the interspecific differences in 
leaf anatomy with the integrated leaf architectural traits such as 
M A and g m ; (iii) to quantify the distribution of overall photosyn- 
thetic limitation among biochemistry, mesophyll diffusion, and 
stomata; and (iv) to quantify the resistance that each anatomical 
component exerts on the diffusion of C0 2 inside the leaf. 



Material and methods 

Plant material 

Fifteen taxa of different growth form and leaf longevity were selected 
for the study to obtain an extensive range of variation in leaf mor- 
phology and anatomy (Supplementary Table SI at JXB online). Five 
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species were annual herbs {Capsicum annuum, Helianthus annuus, 
Phaseolus vulgaris, Spinacea oleracea, Ocimum basilicum) and the 
rest were broad-leaved trees: four deciduous {Acer negundo, Alnus 
subcordata, Betula pubescens, Catalpa speciosa), one semi-deciduous 
{Quercus brant ii) and five evergreens (Quercus ilex, Citrus reticulata, 
Ficus elastica, Pittosporum tobira, Washingtonia filifera). All species 
were dicots, except for the palm Washingtonia filifera. 

All plants were grown either from commercial seed or from 
seeds collected in the field, except for F. elastica where rooted cut- 
tings of a single mother plant were used. The plants were grown in 
a growth room with a 10 h photoperiod, a day/night temperature 
of 24/18 °C, 60% air humidity, and a constant photon flux density 
of 350 umol rrT 2 s~' at plant level provided by Philips HPI-T Plus 
400 W metal halide lamps. The daily integrated incident quantum 
flux density was 12.6mol nr 2 d~'. The growth substrate was a 1/1 
mix of quartz sand and standard potting soil (Biolan Oy, Finland) 
including slow-release NPK (3/1/2 ratio) fertilizer with microele- 
ments, and the plants were irrigated daily to soil field capacity. The 
size of the pots varied between 1 and 5 1 depending on plant age 
and size. In all cases, fully developed young (current-season leaves 
in evergreens) leaves were used for the measurements. In herbs, the 
plants were measured 1 month after seed germination, whilst woody 
species were measured on the second growing year. All physiological 
and structural analyses were replicated with at least three independ- 
ent plants per taxa. 



Foliage gas exchange and fluorescence measurements 

Attached leaves were used for simultaneous leaf gas-exchange 
and chlorophyll-fluorescence measurements using a portable 
gas exchange fluorescence system GFS-3000 (Walz, Effeltrich, 
Germany) equipped with a leaf chamber fluorometer with an 8 cm 2 
cuvette window area. Light was provided by the LED light source 
of the leaf chamber fluorometer (10% blue and 90% red light) and 
the humidity was controlled by a built-in GFS-3000 humidifier. Use 
of a certain fraction of blue light is routinely used in portable pho- 
tosynthesis devices to induce stomatal opening. Although blue light 
is absorbed more strongly by the upper leaf layers and may lead 
to discrepancies among photosynthesis and fluorescence profiles 
(Evans and Vogelmann, 2006), thereby altering g m estimations by 
the combined gas-exchange/fluorescence techniques (Loreto et al., 
2009), the amount of blue light used in our study was small and the 
expected effect minor. 

The standard conditions for leaf stabilization in the cuvette 
were: leaf temperature of 25 °C, saturating quantum flux density 
of 1500 umol nr 2 s~', and C0 2 concentration in the cuvette (C a ) 
of 385 umol C0 2 mol air -1 . Once the steady-state conditions were 
reached, typically 15-20min after clamping the leaf in the cuvette, 
C0 2 response curves of net assimilation (^4 N ) were measured. First, 
C a was lowered stepwise from 385 to 50 umol C0 2 mol air -1 and then 
raised again to 385 umol C0 2 mol air -1 , and the leaf was kept at this 
C a until the original A N value was achieved. Next, C a was increased 
stepwise from 385 to 1500 umol C0 2 mol air 1 and returned again 
to 385 umol C0 2 mol air -1 . In all cases, measurements of ^4 N and 
steady-state fluorescence yield (F s ) were recorded after the gas- 
exchange rates stabilized at the given C a . After recording the A K 
value, a flash of saturating white light was given to determine the 
maximum fluorescence yield in light-adapted state (FJ). After com- 
pletion of the C0 2 response curves, the light was switched off and 
respiration rate in the dark (R d ) was determined. In calculations of 
A-n, R d , and intercellular C0 2 concentration (Q), corrections for 
the diffusion leakage of C0 2 into and out of the leaf chamber were 
included as described by Flexas et al. (2007). 



Measurements of leaf optical properties 

Leaf transmittance and reflectance measurements were conducted 
with a spectrometer (AvaSpec-2048-2; Avantes, Apeldoorn, The 
Netherlands) using an integrating sphere (ISP-80-8-R; Ocean 



Optics, Dunedin, FL, USA). Leaf absorptance (a) was calculated 
as 1 minus the sum of reflectance and transmittance. Three leaves 
of each species were measured, and within each leaf, three rep- 
licate measurements were made. Average absorptance across the 
400-700 nm region was used to characterize the fraction of incident 
photosynthetically active radiation absorbed by the leaf. 



Anatomical measurements 

After the gas-exchange measurements, lxl mm pieces were cut 
between the main veins from the same leaves for anatomical meas- 
urements. Leaf material was quickly fixed under vacuum with 4% 
glutaraldehyde and 2% paraformaldehyde in 0. 1 M phosphate buffer 
(pH 7.2). Afterwards, the samples were fixed in 1% osmium tetroxide 
for 1 h and dehydrated in a graded ethanol series, followed by wash- 
ing three times in propylene oxide. The dehydrated segments were 
embedded in Spurr's resin (Monocomp Instrumentation, Madrid, 
Spain) and cured in an oven at 60 °C for 48 h. Semi-thin (0.8 um) and 
ultrathin (90 nm) cross-sections were cut with an ultramicrotome 
(Reichert & Jung model Ultracut E). Semi-thin cross-sections 
were stained with 1% toluidine blue and viewed under an Olympus 
BX601ight microscope. Photos were taken at 200x and 500x magni- 
fication with a digital camera (U-TV0.5XC; Olympus) to measure 
the leaf thickness and thickness of the palisade and spongy tissue 
layers (Supplementary Fig. S1A-C). Ultrathin cross-sections for 
transmission electron microscopy (TEM H600; Hitachi) were con- 
trasted with uranyl acetate and lead citrate. Photos were taken at 
2000 x magnification (Supplementary Fig. S1D-F) to measure the 
length of mesophyll cells and chloroplasts adjacent to intercellular 
air spaces and chloroplast width and thickness, and the volume frac- 
tion of intercellular air space calculated as: 



where 2S S is the total cross-sectional area of mesophyll cells, IV is 
the width of the section, and ? mes is the mesophyll thickness between 
the two epidermises. Mesophyll (S m /S) and chloroplast (SJS) sur- 
face area exposed to intercellular air spaces per leaf area were cal- 
culated separately for spongy and palisade tissues as described by 
Evans et al. (1994) and Syvertsen et al. (1995). The curvature correc- 
tion factor was measured and calculated for each species according 
to Thain (1983) for palisade and spongy cells by measuring their 
width and height and calculating an average width/height ratio. The 
curvature factor correction ranged from 1.16 to 1.4 for spongy cells 
and from 1.4 to 1.5 for palisade cells. All parameters were analysed 
at least in four different fields of view and at three different sections. 
Weighted averages based on tissue volume fractions were calculated 
for SJS and SJS. 

r cw and cytoplasm thickness (T cyt ) were measured at 20 000^0 
000x magnification depending on the species (Supplementary Fig. 
S1G-I). Three different sections per species and four to six differ- 
ent fields of view were measured for each anatomical characteristic. 
Micrographs were selected randomly in each section and T cvl was 
measured for two to three cells per micrograph. Ten measurements 
for spongy tissue and ten for palisade parenchyma cells were made 
for each anatomical trait, and weighted averages based on tissue vol- 
ume fractions were calculated. All images were analysed with Image 
analysis software (ImageJ; Wayne Rasband/NIH, Bethesda, MD, 
USA). 



M A and leaf density 

The leaves were scanned at 300 dpi, and then oven dried at 70 C for 
48 h and their dry mass was estimated. Leaf area was determined 
from the images with Image J. From these measurements, M A was 
calculated. Using the estimates of leaf thickness from anatomical 
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measurements, leaf density (Z) L ) was calculated as M A per unit leaf 
thickness (Niinemets, 1999). 

Estimation ofg m and model parameters Farquhar et al. (1980) 
by combined gas-exchange/fluorescence method 
The chloroplastic hypothetical C0 2 compensation point (T*) in the 
absence of R<i was calculated from the Rubisco specificity factor 
(S ao ) as: 



parameter significantly improves the model fit. Estimates of J max , 
F cmax , and g m were derived from fitting A N -Q curves as previously 
described. Values of the Michaelis-Menten constant for C0 2 (K c ), 
and oxygen (K a ) and their temperature responses used for these esti- 
mations were from Bernacchi et al. (2002). T* was calculated accord- 
ing to Eqn 2, and R d by gas exchange measurements at 385 umol 
C0 2 mol air -1 . At least three plants per species were used to estimate 
g m . The same leaves were used for estimation of g m by the Ethier and 
Livingston (2004) and Harley et al. (1992) methods. 



E*= 0.5 O/SC/O 



(2) 



using the average values for S a o reported by Galmes et al. (2005) 
for each different leaf habit (Supplementary Table S2 at JXB 
online). A sensitivity analysis showed that the precise value of T* 
within the reported range did not significantly affect the g m estimates 
(Supplementary Table S3A at JXB online). 

From chlorophyll fluorescence measurements, the actual photo- 
chemical efficiency of photosystem II (<t> PSII ) was determined from 
F s and the maximum fluorescence yield during a light-saturating 
pulse of 4500 umol nr 2 s -1 CF m ') following the method of Genty et al. 
(1989): 



(3) 



The linear electron transport rate on the basis of chlorophyll fluo- 
rescence (/ F ) was then calculated as: 



/ F - cDpgjjgae psii 



(4) 



where Q is the photosynthetically active quantum flux density, a 
is the leaf absorptance, and e PSII is the fraction of light absorbed 
by PSII. As routinely assumed, e Psn was taken as 0.5 (Loreto et al., 
1 994; Niinemets et al. , 2005). 

Furthermore, the g m to C0 2 was estimated according to Harley 
etal. (1992) as: 



g m modelled from anatomical characteristics 

The one-dimensional gas diffusion model of Niinemets and 
Reichstein (2003a) as applied by Tosens et al. (2012a) was employed 
to estimate the share of different leaf anatomical characteristics in 
determining g m . g m as a composite conductance for within-leaf gas 
and liquid components is given as: 



1 



H ■ gliq 



(6) 



where g- ias is the gas phase conductance inside the leaf from subs- 
tomatal cavities to outer surface of cell walls, g iiq is the conductance 
in liquid and lipid phases from outer surface of cell walls to chlo- 
roplasts, R is the gas constant (Pa m 3 Kr 1 moL 1 ), T k is the absolute 
temperature (K), and His the Henry's law constant (Pa m 3 mol -1 ). g m 
is defined as a gas-phase conductance, and thus HI(RTy), the dimen- 
sionless form of Henry's law constant, is needed to convert g^ q to 
corresponding gas-phase equivalent conductance (Niinemets and 
Reichstein, 2003a). In the model, the gas-phase conductance (and 
the reciprocal term, r ias ) is determined by average gas-phase thick- 
ness, AL ias , and gas-phase porosity,y5 as (fraction of leaf air space): 



1 



-^a " f\a.\ 



(7) 



in 



r*(/ F + 8(A + * L )) 

J F -4(A N +R L ) 



(5) 



where R L is the respiration rate in the light. In this study, R d was 
used as a proxy for R L (Galle et al., 2009). In other studies, half R d 
has been used (Piel et al., 2002; Niinemets et al., 2005). However, 
as shown in Supplementary Table S3B, no significant differences 
in g m were found when using the proxy for R L , and hence we con- 
cluded that selection of the appropriate value for R L is not a critical 
issue for our g m estimates, confirming a previous sensitivity analysis 
(Niinemets et al, 2006). 

The obtained values of g m were used to transform the A^-Q 
response curves into A^ versus C c response curves as C c = Q - A^/g m . 
Finally, Farquhar et al. (1980) model parameters, the maximum 
velocity of carboxylation (K cmax ) and the capacity for photosyn- 
thetic electron transport (/ max ) on the basis of C c were calculated 
according to Bernacchi et al. (2002). Three replicates estimates of g m 
were available for every species. 

Estimation ofg m from gas exchange measurements only: the 
curve-fitting method 

The curve-fitting method introduced by Ethier and Livingston 
(2004) as applied by Niinemets et al. (2005) was used to obtain an 
alternative estimate of g m . This method is based on changes in the 
curvature of ^4 N versus Q response curves due to finite g m such that 
the Farquhar et al. (1980) model based on Q imperfectly fits the 
data (Ethier and Livingston 2004). Thus, including g m as a fitted 



where is the diffusion path tortuosity (m nr 1 ) and Z) a (m 2 s" 1 ) 
is the diffusion coefficient for G0 2 in the gas phase (1.51 x 10 -5 at 
25 °C). AL ias was taken as half the mesophyll thickness. The partial 
determinants of the liquid-phase diffusion pathway (the reciprocal 
term r b where i stands either for cell wall, cytosol, or stroma con- 
ductance) were calculated as: 



1 _r ty D, 
i] 



A Li 



(8) 



where ALj (m) is the diffusion path length in the correspond- 
ing component of the diffusion pathway, p { (m 3 m~ 3 ) is its effective 
porosity, and D w is the aqueous phase diffusion coefficient for C0 2 
(1.79 x 10~ 9 m 2 s -1 at 25 °C). The dimensionless factor r f i accounts 
for the reduction of Z) w compared with free diffusion in water, and 
was taken as 1.0 for cell walls (Rondeau-Mouro et al., 2008) and 0.3 
for cytosol and stroma (Niinemets and Reichstein, 2003b). In addi- 
tion, f fi values for cytosol and stroma were estimated using a least- 
squares iterative analysis to assess the sensitivity of g m to values of 
r si (Supplementary Figs S2 and S3 at JXB online). In this analysis, 
r fi was allowed to vary between 1 and 0.05, and the values of r (i were 
varied within this range to minimize the difference between meas- 
ured and modelled g m . Whilst this approach improved the agreement 
between modelled and measured g m , the extreme values obtained for 
r si seemed unrealistic (Supplementary Figs S2 and S3). p { was set 
to 1.0 for cytosol and stroma. There are no direct measurements 
of cell wall porosity, but it has been suggested that this parameter 
might vary with T cw among species (Terashima et al., 2006; Evans 
et al., 2009; Tosens et al., 2012b). Therefore, given the heterogeneous 
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series of species used in this experiment, p { was estimated using a 
least-squares iterative analysis assuming a hypothetical relation- 
ship between porosity and T cw as described by Tosens et al. (2012b). 
Again, a least-squares iterative approach was employed to get the 
best fit between measured and modelled g m . The Pi range in this 
analysis was fixed at 0.028 (Tosens et al., 2012b) for the thickest to 
0.3 (Nobel, 1991) for the thinnest cell walls (Supplementary Table 
S5 at JXB online). We used an estimate of 0.0035 m s -1 for both 
plasma membrane conductance (g pi ) and chloroplast envelope con- 
ductance (gem) as previously suggested (Evans et al., 1994; Tosens 
et al., 2012a). 

Carbonic anhydrase in cytosol and chloroplasts could facilitate 
the diffusion of C0 2 through the liquid phase. However, there is lit- 
tle evidence for the involvement of carbonic anhydrase in g m and A K 
(reviewed by Flexas et al., 2008, 2012). Therefore, following Tosens 
et al. (2012a), we did not include the potential effect of carbonic 
anhydrase in our analysis. 

In previous studies, we scaled the total liquid-phase diffusion 
conductance by SJS ratio (Tosens et al., 2012a) that determines the 
number of parallel diffusion pathways from outer surfaces of cell 
walls to chloroplasts. 



respectively, both calculated by Eqn 8. For r cyt u the diffusion path- 
way length, AL cyt>1 , is given as the average distance between the chlo- 
roplasts and cell wall in cell wall areas lined by chloroplasts (Fig. 1), 
whilst for r st j, AL { , was taken as half of the chloroplast thickness, 
AT C J2. For the cell wall portions without chloroplasts, the partial 
conductance, g cA2 , is given analogously as: 



gcel,2 



r cyi,2 + r env + ?st,2 



(ID 



where r cyU2 is the cytosolic resistance from interchloroplastic cell 
wall portions towards the chloroplast and r s , 2 is the stromal conduct- 
ance in a direction parallel with the cell wall (Fig. 1). The diffusion 
path length for j' cyt>2 (Eqn ), AL cyt2 , is driven both by the distance 
between the neighbouring chloroplasts, chloroplast thickness, and 
chloroplast distance from the cell wall and was approximated as: 



AL, 



tyt,2 " 



AL, 



-AL, 



'cyl.l 



A La 



(12) 



giiq = - 



r cw + 'pi + 'cyt + r en + 'st ) S 



(9) 



Although, cell wall and plasmalemma resistances actually scale 
with the SJS ratio, use of SJS has been deemed to be appropri- 
ate, as SJS is generally close to the SJS ratio (Scafaro et al., 201 1; 
Peguero-Pina et al., 2012), i.e. there is little cell wall area free of chlo- 
roplasts. Even if SJS is significantly smaller than SJS, the cytosolic 
distance between the neighbouring chloroplasts is generally large 
and this can still constrain the diffusion flux in interchloroplastial 
areas of cell wall (locally large cytosol conductance, g cyt ; Fig. 1). 
However, the significance of the ;- cyt depends on the other parts of 
the diffusion pathway as well. 

To explicitly assess the importance of diffusion of C0 2 through 
interchloroplastial areas, we considered two different pathways of 
C0 2 inside the cell, one for cell wall parts with chloroplasts and 
the other for interchloroplastial areas as described by Tholen et al. 
(2012). For exposed cell wall portions lined with chloroplasts, the 
partial liquid phase conductance, g^u inside the cell is given as: 



1 



(10) 



cyt.l 



'st.l 



where 



cyt,l 



and r st j are cytosolic resistance from the plasma- 



lemma inner surface to the outer surface of chloroplasts and the 
stromal resistance in the direction perpendicular to cell wall (Fig. 1), 



where AL chl is the distance between the neighbouring chloro- 
plasts. AZ cyt>2 was calculated as a harmonic average, which more 
correctly represents the diffusion pathway of r cyt2 . Finally, the dif- 
fusion pathway length for r sl2 was taken as a quarter of the chlo- 
roplast length. 

Considering further that the fraction of exposed cell wall area 
lined with chloroplasts is given by SJS m and the fraction free of 
chloroplasts as 1 - SJS m , the total cellular conductance (sum of par- 
allel conductances) is given as: 



> eel, 2 



(13) 



Total liquid phase conductance from the outer surface of cell 
walls to carboxylation sites in the chloroplasts is the sum of serial 
conductances in the cell wall, plasmalemma, and inside the cell: 



gliq : 



(14) 



'™ + r„, + r 



cel,tot j 



Alternatively, the total cellular diffusion pathway can be consid- 
ered to consist of two parallel pathways from the outer surface of 
the cell walls to the chloroplasts, one pathway corresponding to the 
diffusion flux through cell wall areas lines with chloroplasts and the 
other without chloroplasts: 



CM 



Path 1 



Stroma 
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is 

Q. 



Path 1 
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Cell wall 



Chloroplast 
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ALcytl 



AL 
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Fig. 1. Illustration of the diffusion pathway in exposed cell wall areas lined with chloroplasts (path 1) and interchloroplastial areas (path 
2). The diffusion pathway in leaf lipid and liquid phases includes cell wall, plasmalemma, cytosol (shown by red arrows), chloroplast 
envelope membranes, and chloroplast stroma (shown by dark green arrows). The effective diffusion path length in cytosol along path 
1 is taken as the average distance of chloroplasts from the cell wall, AL cyt1 , whilst the diffusion pathway length in interplastidial areas is 
determined by the distance between the chloroplasts and A/_ cyt1 (Eqn 12). 
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gliq : 



i) S (j- cw +;- pl + j- celi2 ) S 



(15) 



were revealed by Tukey analyses (P <0.05). These analyses were per- 
formed with the IBM SPSS statistics 19.0 software package (SPSS). 



Although Eqns 14 and 15 are conceptually different, the values of 
conductances calculated by both equations were very similar, differ- 
ing at most by 4%. In the current study, we have used Eqn 14. 



Analysis of quantitative limitations on A N 

To separate the relative controls on A N resulting from limited 
stomatal conductance (/ s ), mesophyll diffusion (/ m ), and limited 
biochemical capacity (/ b ) (4+4n + 4> =1 )> we used the quantitative 
limitation analysis of Jones (1985) and implemented by Grassi and 
Magnani, (2005). The limitations of the different components were 
calculated as: 



I, 



gtot/gs 'dA^/dC c 
g m + dA N I dC c 



. gtotlgm'dAnldCc and 
gtot + dAx I dC c 



k = - 



gtot 



gtot 



■dA*ildC r 



(16) 



(17) 



(18) 



where g s is the stomatal conductance to C0 2 , g m was according 
to Harley et al. (1992, Eqn 5), and g, ot is the total conductance to 
C0 2 from ambient air to chloroplasts (sum of the inverse serial con- 
ductances g s and g m ). 8A^/8C C was calculated as the slope of A^-C c 
response curves over a C c range of 50-100 umol mof 4 . At least three 
curves per species were used, and average estimates of the limita- 
tions were calculated. 



Quantitative analysis of partial limitations ofg m 

The determinants of g m were divided between the component parts 
of the diffusion pathway (Eqns 6-8). The proportion of g m deter- 
mined by limited gas-phase conductance (/ ias ) was calculated as: 



(19) 



The share of g m by different components of the cellular phase 
conductances (/;) was determined as: 



' S 



(20) 



where /; is the component limitation in the cell walls, the plasma- 
lemma, and inside the cells, and g ; refers to the component diffusion 
conductances of the corresponding diffusion pathways. To determine 
the limitations derived from the different components inside the cell 
(cytoplasm, chloroplast envelope, and stroma), weighted limitations 
of both pathways, the fraction of exposed cell wall area lined with 
chloroplasts and the fraction free of chloroplasts, were used. 



Statistical analyses 

Regression and correlation analyses were conducted using the Sigma 
Plot 10.0 software package (SPSS; Chicago, IL, USA). Univariate 
analysis of variance was performed to reveal differences between 
species in the studied characteristics. Differences between means 



Results 

Leaf structural and anatomical traits 

M A varied sixfold (20-123 g rrT 2 ) (Supplementary Table S4 at 
JXB online). The variation in leaf thickness was 3.7-fold with 
Acer negundo having the thinnest (123 um) and F. elastica the 
thickest (459 um) leaves. Spongy mesophyll thickness varied 5.2- 
fold, and palisade mesophyll thickness 2.5-fold (Supplementary 
Table S4). Generally, the palisade tissue comprised approxi- 
mately 40%, and spongy tissue approximately 60% of total 
mesophyll, except for some species as F. elastica with 75% and 
W. filifera with 100% of spongy tissue. The variation in D L was 
6.4-fold with Phaseolus vulgaris having the least dense (0. 1 1 g 
cnT 3 ) and Q. ilex the most dense (0.70 g cnT 3 ) leaves. M A exhib- 
ited a significant positive correlation with D L (Supplementary 
Fig. S4 at JXB online), but was weakly correlated with leaf 
thickness (^=0.27, P <0.05; data not shown). Therefore, the 
variation in M A was mainly attributed to the leaf density. 

Among the leaf ultrastructural characteristics estimated 
from transmission electron micrographs (Supplementary 
Tables S4 and S5, and Supplementary Fig. S1D-I), SJS var- 
ied 3.3-fold across all species (14.4-40 m 2 nr 2 ) and SJS varied 
2.7-fold (6-19.7 m 2 nT 2 ). SJS m varied between 0.31 (Citrus 
reticulata) and 0.74 (O. basilicum). For T CVI (Supplementary 
Fig. S1G-I), 4.8-fold variation was observed between all 
species (113. 6-543. 7nm). Herbaceous species exhibited the 
thinnest cell walls together with Catalpa speciosa, whilst ever- 
greens had the thickest cell walls with the maximum value of 
543.7 nm observed in Pittosporum tobira. 

Estimation ofg m with different methods 

The values of g m calculated according to the methods of 
Harley et al. (1992) and Ethier and Livingston (2004) were 
strongly correlated (Supplementary Fig. S5 at JXB online, 
r 2 =0.80). However, the Harley et al. -based estimates exhibited 
the smallest average coefficient of variation for independent 
estimates within a species and therefore we report the data 
obtained with this method only. 

Mesophyll conductance calculated by the method of 
Harley et al. (1992) varied 24-fold across all species. H. ann- 
uus showed the maximum values and Citrus reticulata the 
minimum values of g m . The minimum value for the coefficient 
of variation in g m was 1.9% (Pittosporum tobira), whilst the 
maximum value was 32.9% (Q. ilex). The average of the coef- 
ficient of variation for all species was 16.5%. 

g m in relation to physiological characteristics 

Net assimilation rate correlated positively with g s and g m 
(Supplementary Fig. S6 at JXB online). Q-C c ranged from 
240 to 112 umol mok 1 in woody deciduous and evergreen 
species, and had lower values (40-67 umol moF 1 ) in herbs. 
C r C c decreased with increasing g m (Supplementary Fig. S7 
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at JXB online). This relationship was qualitatively identical 
when g m was expressed on the leaf area or dry mass basis 
(data not shown). 

g m in relation to leaf structural and anatomical traits 

g m per dry mass was negatively associated with M A (r 2 - 0.85, P 
<0.0005; data not shown). g m per unit leaf area or per unit dry 
mass (data not shown) was not correlated with SJS, reflecting 
the circumstance that SJS was almost invariable, between 16 
and 24 m 2 nr 2 across the species. SJS was not significantly cor- 
related with g m (Fig. 2A, P >0.13). However, a positive correla- 
tion between g m and SJS was observed when the species with the 
largest T cw (Pittosporum tobira and Q. brantii) were not included 
in the correlation (t 2 -0J7, P <0.0001; data not shown). 

The positive and significant correlation (r 2 =0.84, P <0.001) 
between g m and (SJS)/D L suggested the importance of the 
anatomical components to the internal diffusion of C0 2 
(Fig. 2B). Moreover, the negative and significant relationship 
observed between gJ(SJS) and r cw showed the importance 
of r cw in affecting g m (Fig. 2C). 

g m calculated from anatomical variables 

Using the leaf anatomical traits measured, g m was modelled 
and compared with g m measured by the method of Harley 
et ol. (1992). A good positive linear relationship between mod- 
elled and measured g m was observed (^=0.90, P <0.0001; 
Fig. 3). However, the slope was different from unity, so that the 
g m modelled tended to be overestimated in species with high 
M A and underestimated in species with low M A . g m values cal- 
culated by the model based on leaf anatomy ranged between 
0.217 and 0.056mol nr 2 s" 1 . H. annum showed the largest and 
W. filifera the smallest values of g m . The coefficient of intraspe- 
cific variation in g m estimates for different replicates was lower 
than for the experimental estimations, being between 1.2% 
(Betula pubescens) and 22% (Pittosporum tobira). 



Overall importance of g m 

According to quantitative limitations analysis of A^, sto- 
matal openness and g m restricted the photosynthetic capac- 
ity to a similar percentage, 19-65% and 13-64%, respectively. 
However, the biochemical limitations were lower than the 
stomatal and mesophyll limitations, being between 6 and 33% 
(Fig. 4A-C). Both the stomatal and biochemical components 
tended to be more important in species with non-sclerophytic 
leaves (low M A ), whilst mesophyll diffusion limitation was most 
significant in species with high M A (Fig. 4). Thus, herbaceous 
plants showed the maximum values for stomatal limitations, 
whilst the maximum mesophyll limitations were observed in 
evergreen species with more robust foliage structure. 

Limitation of g m due to individual components of the 
diffusion pathway 

From the different components of the whole diffusion path- 
way of C0 2 , the percentage limitations of g m were estimated 
(Fig. 4D-I). Intercellular air spaces represented a smaller resist- 
ance to the C0 2 diffusion (4-22%) than the cellular phase, 
because the rate of C0 2 diffusion in air was larger than in 
water. In the cellular phase, the cell walls appeared to be the 
most important factor that limited the internal diffusion of C0 2 
in the species that presented a high M A . However, the plants 
with low M A that presented a low percentage of limitation of 
g m by the cell wall revealed a higher limitation by the stroma of 
around 43%. On the other hand, the plasmalemma and chloro- 
plast envelope accounted for only up to 8% of the limitation. 

Discussion 

Values of g m in a range of species exhibiting different 
foliage morphologies 

The range of g m values observed in our study is representa- 
tive of the whole range of g m values described so far in large 
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Fig. 2. Correlations of mesophyll diffusion conductance (g m ) determined according to Harley et al. (1 992) with the surface area of 
chloroplasts exposed to intercellular airspaces per unit leaf area (S/S) (A), mesophyll diffusion conductance with chloroplast surface 
area per leaf density ((S C /S)/D|J (B), and mesophyll diffusion conductance per S c /S (gJ{S c /S)) with T" cw (C). In the main panels, the data 
were fitted by linear (B) and non-linear (C) regressions in the form y=ae~bx. In the inset, the data were fitted by linear regression. Different 
species are represented as: herbs (circles), woody deciduous and semi-deciduous species (triangles), and woody evergreen species 
(squares). Values are means ±standard error (SE) of three to four replicates per species. 
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g modeled (mol CO, m"V) 

Fig. 3. The relationship between mesophyll diffusion conductance 
(gj measured with Harley ef al. method and g m modelled with 
anatomical parameters (Eqn 6-1 5). Values are means +SE of three 
replicates per species. Symbols are the same as in Fig. 2. The 
data were fitted by linear regression. Broken lines correspond to 
the 1 :1 relationship. 

literature-based datasets, except that the maximum g m values 
found in the present study were somewhat lower than reported 
previously (Flexas et al, 2008; Warren, 2008, Niinemets et al, 
2009a). These relatively low maximum values were possibly 
due to moderate growth light intensity compared with full 
sun (Piel et al., 2002; Niinemets et al., 2009a). This explana- 
tion is consistent with the observations of a significant num- 
ber of chloroplasts not closely facing the cell walls (Fig. 2A) 
and relatively low ratios of chloroplast exposed to mesophyll 
exposed cell wall surfaces {SJS m ) (Supplementary Table S4), 
both being traits that depend on the growth light environ- 
ment (Terashima et al, 2006). 

Relationship ofg m to leaf anatomy and its importance 
in limiting photosynthesis 

As in previous studies, g m showed a high degree of correla- 
tion with several leaf anatomic characteristics, notably a nega- 
tive correlation with M A (Flexas et al., 2008; Niinemets et al., 
2009a,6) and a positive correlation with SJS (Evans et al., 
1994, 2009). The M A effect on g m supports the idea that g m 
depends on species differences in leaf density, as density was 
positively correlated with M A , whilst leaf thickness showed a 
weak correlation. The photosynthetic capacity was also sig- 
nificantly and positively correlated with g m as demonstrated 
previously (reviewed by Flexas et al, 2008; Niinemets et al., 
2009a). Overall, these results suggest that, in species with high 
M A , photosynthesis is more limited by g m , as indirectly sup- 
ported by the negative effect of leaf density on g m (Niinemets, 
1999) and more directly evidenced by the fact that they present 
higher values of C r C c (Warren, 2008; Niinemets et al, 2009a). 



The relative contribution of g s , g m , and photosynthetic 
biochemistry to total photosynthesis limitation (follow- 
ing Grassi and Magnani, 2005) was variable and depended 
on leaf structural characteristics, i.e. M A (Fig.4A-C). At a 
typical operating C0 2 concentration, the biochemical limi- 
tations of photosynthesis decreased from a maximum of 
approximately 33% at low M A to minimum values as M A 
increased, whilst, in parallel, mesophyll diffusion limita- 
tions increased from a minimum of approximately 15% to 
maximum values up to 65%. Stomatal limitations showed 
a less clear variation with M A . Overall, these data demon- 
strated that species with low M A showed a notable coordi- 
nation of the limiting factors for photosynthesis, i.e. they 
were similarly co-limited by stomatal, mesophyll, and bio- 
chemical limitations. In contrast, species with high M A were 
mostly limited by mesophyll (on average by 57%) and sto- 
matal (30%) diffusion, and were less limited by biochemis- 
try (13%). This is consistent with the idea that species with 
thicker and denser leaves, e.g. evergreen trees, are more lim- 
ited by g m than species with thinner leaves (Galmes et al, 
2007; Niinemets et al, 201 1). 

Key structural factors regulating differences in g m 
between distant leaf structures 

The fact that g m and mesophyll diffusion limitations were 
strongly correlated with M A suggested that interspecific varia- 
tions in g m are driven by leaf structural characteristics. Among 
the key structural traits suggested to limit C0 2 diffusion the 
most are the traits that alter effective diffusion path length 
and area for diffusion, in particular T cw , and chloroplast dis- 
tribution along the exposed mesophyll cell wall (Fig. 2; Evans 
et al, 2009), although the role of other variables, such as leaf 
porosity, and the path lengths for C0 2 through the plasma- 
lemma and chloroplast envelope membranes, cytosol, and 
stroma cannot be ruled out (Evans et al, 2009). In the present 
study, we modelled g m considering all major leaf structural 
traits as described by Tosens et al. (2012a). A high signifi- 
cant positive correlation between measured and modelled g m 
estimates was found (Fig. 3). This correlation supports the 
view that at least a significant proportion of the interspecific 
variations in g m is somehow related to differences in the thick- 
ness of the structures involved in C0 2 diffusion, as well as to 
the number of parallel C0 2 diffusion pathways determined 
by SJS. 

Despite the high correlation, the slope of the relationship 
was not unity, so that the biggest discrepancies between meas- 
ured and modelled estimates of g m were found at the higher 
and lower ends of g m . A similar discrepancy was observed in 
different Australian sclerophyll species occurring in the field 
under different soil nutrients and water availabilities, espe- 
cially at high values of g m (Tosens et al, 2012b). 

This strong discrepancy between measured and mod- 
elled values may arise from the inherent uncertainties asso- 
ciated with both estimates. As for the Harley et al. (1992) 
approach, besides the small variability in the estimates asso- 
ciated with uncertainties in the exact values of R L and T* 
(Supplementary Table S3), it has recently been shown that g m 
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Fig. 4. Quantitative limitation analysis of photosynthetic C0 2 assimilation and mesophyll conductance to C0 2 (g^) in relation to the leaf 
dry mass per area {M^. The stomatal (A), mesophyll (B), and biochemical (C) limitations of photosynthetic assimilation were calculated 
according to Eqns 1 6-1 8. Limitations analyses were based on a chloroplastic C0 2 concentration (C c ) range of 50-1 00 (imol mor 1 . 
Quantitative limitations of g m due to different anatomical components of the diffusion pathway were calculated using leaf anatomical 
characteristics (Eqns 19 and 20). The relative C0 2 diffusion limitations separated were: intercellular spaces (D), cell wall thickness (E), 
cytoplasm (F), plasmalemma (G), chloroplast envelope (H), and chloroplast stroma (I). 



cannot be considered as a purely diffusional component, but 
instead intrinsically includes a flux-weighted quantity related 
to the amount of respiratory and photorespiratory C0 2 from 
the mitochondria diffusing towards the chloroplasts (Tholen 
et al, 2012). Concerning the anatomically based model used 
here, the precise outputs largely depend on a number of varia- 
bles assumed as constants or inferred indirectly. For instance, 
the reduction in D w compared with free diffusion in water 
(r f ;) was considered constant for all species, although differ- 
ent for cell wall and intercellular components. Both g pl and 
g env were also taken as constant, whilst cell wall porosity (p L ) 
was indirectly estimated from r cw using an empirical equa- 
tion. There is not sufficient knowledge for the actual values 
of all these parameters, and they may vary among species, 
hence contributing to most of the observed slope discrep- 
ancy. It can be seen, for instance, that the difference between 
measured and modelled g m almost disappeared when the r fi 
values were calculated using a least-squares iterative analysis 



(Supplementary Fig. S2). However, this 'perfect correspond- 
ence' is bound to some probably non-realistic r {l values as 
low as 0.05. Moreover, values of g m have been modelled con- 
sidering C0 2 diffusivities in the different media involved — 
assumed to be either 'pure' air, lipid, or aqueous phases with 
fixed thicknesses, whilst, in most cases, determination of the 
thickness of the given phase is not that straightforward. Also, 
we assumed no facilitation mechanism that could improve 
the diffusivities in lipid and aqueous phases. Among these, 
membrane-bound aquaporins (Uehlein et al., 2003, 2008; 
Hanba et al, 2004; Flexas et al, 2006) and cytosol and stro- 
mal forms of carbonic anhydrases (Price et al, 1994; Gillon 
and Yakir, 2000) are likely candidates (Terashima et al., 
201 1). For instance, allowing g pl and/or g env to vary within the 
range of published values (Evans et al, 2009) also results in a 
better agreement between the measured and modelled values 
(Supplementary Fig. S8 at JXB online). In summary, current 
uncertainties about the actual values of these parameters and 
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their variability among species preclude the development of 
a truly predictive anatomically based model for g m . However, 
the good correlation, despite the divergent slope, can be taken 
as strong evidence that a substantial part of g m is indeed 
dependent on a series of leaf anatomical features. 

Despite the discussed limitations of the model approach 
used here, the results suggest that chloroplast distribution 
and r cw are the most influential leaf structural characteris- 
tics in setting the limits for g m (Evans et al., 2009; Terashima 
et al., 2011). In particular, a significant positive correlation 
was found between g m and SJS only when species with very 
large T CVi (Q. brantii and Pittosporum tobira) were excluded, 
highlighting the fact that the impact of chloroplast distribu- 
tion on g m became less important as T CVI increased, in agree- 
ment with past suggestions (Terashima et al., 2006, 201 1). 

In addition, a highly significant negative relationship 
was observed between the ratio gJ(SJS) and r cw consider- 
ing all species, similar to that obtained by Terashima et al. 
(2011) pooling literature data. Using a limitation analysis to 
separate the contributions of the components of g m (Eqns 
13 and 14) revealed that, globally, the limitation imposed by 
T cv/ spanned the most, ranging from approximately 4 to 70% 
(Fig. 4E). This was followed by chloroplast stroma, which 
ranged from 4 to 46% (Fig. 4F, I). However, the limitations 
inside the cell (cytosol and stroma) could be underestimated, 
especially in species with high M A , as g m was modelled assum- 
ing that cytosolic and stromal viscosity (r fi ) was constant in 
all species. 

The limitations imposed by intercellular air spaces, the 
plasmalemma, and the chloroplast envelope were much 
smaller than the rest of the diffusion pathway components as 
was observed by Tosens et al. (2012a,6). The fact that the lat- 
ter two components had only a moderate effect on limiting g m 
is in conflict with the observed larger g m changes observed in 
aquaporin mutant plants without any appreciable differences 
in SJS or any other leaf structural characteristic (Flexas 
et al, 2006). This could be due to the fact that the assumed 
values for g pl and g em are constant among species, which may 
not necessarily be the case. Differences of up to four orders 
of magnitude have been reported for C0 2 permeabilities of 
biological membranes. For instance, if the permeability for 
a given species was 0.00002 m s -1 , as found for chloroplast 
envelopes by Uehlein et al. (2008), instead of the 0.0035 m s" 1 
used in the present simulation, the combined limitation to g m 
imposed could be larger than 40% (data not shown). At the 
other extreme, if values were closer to the 0.016 m s _1 reported 
by Missner et al. (2008) for lipid bilayers, the maximum mod- 
elled g m values will be closer to estimates based by the Harley 
et al. (1992) approach (data not shown). Clearly, improved 
knowledge on the actual permeability to C0 2 of biological 
membranes is required to fully understand the basis for the 
regulation of g m . 

Despite these general tendencies, the impact of each spe- 
cific leaf component on g m changed with M A . Specifically, the 
limitations imposed by T cw strongly increased with increasing 
M A , whilst the limitations associated with all the other com- 
ponents decreased with increasing M A . Thus, in species with 
low M A , like annual herbs, about 60% of the total limitation 



to g m is imposed by cytoplasm and stroma, whilst another 
12% is accounted for by the plasmalemma and chloroplast 
envelope. Moreover, in species with thinner leaves, the frac- 
tion of exposed cell wall lined with chloroplasts (g ce i,i) was 
higher, whilst limitations inside the cell through intercho- 
roplastial areas (g ce i, 2 ) were more important in species with 
higher M A (Fig. S9 at JXB online). This suggests that it is in 
such species where facilitating mechanisms (aquaporins, car- 
bonic anhydrases, chloroplast movements, and others) have 
the strongest influence on g m . In contrast, in species with high 
M A , like evergreen sclerophylls, g m is mostly (up to 70%) lim- 
ited by r cw , which is likely to be less variable in the short term, 
and may explain the low photosynthetic capacity displayed 
by these plants even under non-limiting conditions. Possible 
interspecific variation in the role of aquaporins in limiting g m 
is clearly a topic that deserves high priority in future studies. 

In conclusion, the present study showed that mesophyll 
limitations are crucial in determining the maximum photo- 
synthetic capacity when a large range of leaf types are ana- 
lysed collectively. These limitations are variable depending on 
the leaf structural properties, i.e. M A and associated structural 
traits such as leaf density. The variability in mesophyll diffu- 
sion limitations was explained mainly by variations in the rate 
of C0 2 diffusion pathways through cell walls, as well as the 
area for diffusion determined by the chloroplast distribution. 
However, the impact of each component of the diffusion path- 
way largely depended on M A , so that C0 2 diffusion in species 
with thin leaves (e.g. herbs) depends more on membranes and 
aqueous compartments — and is probably more influenced by 
aquaporins and carbonic anhydrases. In contrast, diffusion in 
species with thick leaves is almost fully determined by cell wall 
conductance. Altogether, the variability in g m with M A helps 
explain the worldwide leaf economics spectrum showing a neg- 
ative dependency between photosynthetic capacity and M A . 



Supplementary data 

Supplementary data can be found at JXB online. 

Supplementary Table SI. List of studied species, species 
origin, life form, and leaf longevity. 

Supplementary Table S2. Physiological characteristics 
measured in all studied species. 

Supplementary Table S3. Sensitivity analysis of the influ- 
ence of uncertainties in chloroplastic hypothetical C0 2 com- 
pensation point (T*) and day respiration on the estimation of 
mesophyll conductance (g m ). 

Supplementary Table S4. Leaf dry mass per unit area (M A ), 
leaf thickness (T L ), leaf density (D L ), thickness of mesophyll 
layers, number of palisade cell layers, mesophyll surface area 
exposed to intercellular airspace (SJS), chloroplast surface 
area exposed to intercellular airspace (SJS), and the ratio 
SJS m in all studied species. 

Supplementary Table S5. Cell wall thickness (T cw ), cyto- 
plasm thickness (T cyt ), chloroplasts length (L chi ), chloroplasts 
thickness (T M ), and effective porosity of the cell wall (p { ). 

Supplementary Fig. SI. Representative light micrographs 
at 200 x magnification for Phaseolus vulgaris, Ficus elastica, 
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and Washingtonia filifera, and representative transmission 
electron micrographs at 2000x magnification for Helianthus 
annuus, Acer negundo, and Washingtonia filifera and at 20 
000 x magnification for H. annuus, Alnus subcordata and 
Pittosporum tobira. 

Supplementary Fig. S2. The relationship between meso- 
phyll diffusion conductance (g m ) measured with the Harley 
Harley et al. (1992) method and g m modelled with anatomi- 
cal parameters using least-squares iterative analysis for the r {l 
parameter. 

Supplementary Fig. S3. Effects of the parameter r fi of the 
cytosol and chloroplast stroma on g m modelled from anatom- 
ical characteristics. 

Supplementary Fig. S4. Correlation between leaf density 
(D L ) and leaf dry mass per unit area (M A ). 

Supplementary Fig. S5. Relationship between g m measured 
according to Harley et al. (1992) method versus the Ethier 
and Livingston (2004) method. 

Supplementary Fig. S6. Net photosynthesis rate C4 N ) in 
relation to stomatal (g s ) and mesophyll (g m ) conductance. 

Supplementary Fig. S7. The relationship between g m and 
C0 2 drawdown (Q-Q). 

Supplementary Fig. S8. The relationship between meso- 
phyll diffusion conductance (g m ) measured with the Harley 
et al. (1992) method and g m modelled with anatomical 
parameters using different values for the membrane perme- 
abilities of plasmalemma (g pl ) and chloroplast membrane 
(gem) conductances. 

Supplementary Fig. S9. Quantitative limitation analysis of 
conductance to C0 2 inside the cell (g ce i,tot) calculated on the 
basis of leaf anatomical characteristics. 
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